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The magnetism and electronic structure of Li-doped SnC>2 are investigated using first-principles 
LDA calculations. We find that Li induces magnetism in SnC>2 when doped at the Sn site but 
becomes non-magnetic when doped at the O and interstitial sites. The calculated formation energies 
show that Li prefers the Sn site, in agreement with previous experimental works. The interaction 
of Li with native defects (Sn Vs n and O Vo vacancies) is also studied, and we find that Li not 
only behaves as a spin polarizer, but also a vacancy stabilizer, i.e. Li significantly reduces the defect 
formation energies of the native defects and helps the stabilization of magnetic oxygen vacancies. The 
electronic densities of states reveal that those systems, where the Fermi levels touch the conduction 
(valence) band, are non-magnetic (magnetic). 

PACS numbers: 71.22. +i, 71.55.-i, 75.75.Pp, 61.72.jd, 61.72.jj 



hand Chang et ali22 observed FM in very thin pristine 
Sn02 films which was induced by the presence of oxy- 
gen vacancies located near the film surface. There is also 
a belief that the observed FM in undoped S11O2 mainly 
originates from the bulk double oxygen vacancies^ In 
contrast Wang et alM- studied nanosheets of Sn02 and 
they found that the saturation magnetization of all their 
annealed samples did not feature mono- dependence on 
oxygen vacancies, whereas an Sn vacancy related origin 
was accounted for variations in the magnetization of their 
studied samples. 

Light elements can induce magnetism, but the observed 
magnetism has also a linkage with native defects. How- 
ever, whether magnetism is induced by cation or anion 
vacancies in the pristine host material or doped with 
light elements, the defect formation energies of native de- 
fects, which are important for magnetism, are very high. 
Oxygen vacancies rather than tin vacancies are consid- 
ered as a source of magnetism due to its lower formation 



I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) have recently 
been a major focus of magnetic semiconductor research. 
In the recent past, ferromagnetism (FM) above room 
temperature have been extensively investigated in transi- 
tion metal (TM)-doped wide band gap oxide semiconduc- 
tors, e.g. ZnO 1 ^ and T1O2 — ~— . Among these wide band 
gap semiconductors Sn02 has long been a material of in- 
terest for applied and pure research purposes. In 2002, 
Ogale et al& not only reported room temperature FM in 
Co-doped Sn02, but also very large magnetic moments. 
Later on density functional theory (DFT) calculations 
showed that Sn vacancies, Vs n , can induce magnetism in 
Sn02 and he observed large magnetic moment was at- 
tributed to its influenced Recently, several experimental 
and theoretical reports showed that cation vacancies can 
induce magnetism not only in Sn02, but also in ZnO, 
Ti0 2 , Zr0 2 , In 2 3 ,Ce0 2 HfOA— • 

Following the prediction of magnetism without TM im- 
purities, much interest was also diverted to magnetism in- 
duced by light element doping in host semiconductor ma- 
trices, e.g., K, N, Mg, and C-doped Sn0 2 — ~— ■ There are 
experimental evidences which clearly demonstrate room 
temperature FM induced by light elements ^* -" However, 
the exact nature of magnetism in these semiconductor 
oxides is still under debate. Recent experimental reports 
claim that pristine S11O2 nanocrystalline thin films ex- 
hibit room temperature FM, but the magnetic moment is 
suppressed when doped with Gd^i The presence of large 
amounts of singly ionized oxygen vacancies (Vq), rather 
than Vs n or Vo were found to be responsible for the ob- 
served FM in pristine Sn0 2 thin films^i On the other 



energy^ but neutral Vo does not induce magnetism in 
oxidesi 7 ' 20 i 26 The major issue in magnetism induced by 
vacancies is therefore that these magnetic systems have 
high formation energies 27 i 28 . Consequently, to realize 
defects-driven magnetism experimentally it is essential to 
reduce the defect formation energy of the host material 
(Sn0 2 in our case). We choose Li as a dopant^ that can 
modify the defect formation energies of native defects in 
Sn02- We show that Li can indeed significantly reduce 
the formation energies of various types of defects. This is, 
as far as we know, the first theoretical study performed 
on the interaction of Li with native defects in Sn02- 
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FIG. 1: (color online) A 2 x 2 x 2 supercell of SnC>2. In 
(a) V represents the Sn vacancy, which is fixed at the center 
of the supercell. The positions of Li doped at different Sn 
sites are marked as 1-6, whereas the positions of Li doped at 
O sites are shown by I, II, and III. In model (b) the central 
atom represents the Li atom doped at the Sn site and the O 
vacancies are marked as 1-4. In model (c) Li is doped at the 
O site, marked as Li, and the oxygen vacancies are created at 
the oxygen sites marked as 1, 2, and 3. Big (green) and small 
(blue) balls represent Sn and O atoms, respectively. 



II. COMPUTATIONAL MODELS AND 
METHOD 

S11O2 is known to crystallize in the rutile structure with 
space group P^/mnm or Df% (SG136) under ambient 
conditions^ We used our previously optimized lattice 
parameters.— To simulate the defects a 2 x 2 x 2 supercell 
of SnC>2 containing 48 atoms was employed. We con- 
sidered the following models to investigate the detailed 
energetics and magnetic properties of different defects in 
Sn0 2 . 

a) Li-doped Sn02: Li was doped at the Sn, O, and 
interstitial sites, denoted as Lis n , Lio, and Lij nt , 
respectively. 

b) Vsn+ Lisnj :In this case, Vs n was fixed at the center 
of the supercell, and Li was doped at different Sn j 
sites, with j = 1 — 6 (see Fig.QJa)). 

c) Vsn+ Lioj: In this case, Vs n was fixed at the center 
of the supercell, and Li was doped at different O j 
sites, with j=I, II, III (see Fig. QIa)). 

d) Vsn+ Lij nt : Here Vs n was fixed at the center of the 
supercell and Li was doped at interstitial sites. 

e) V0.7+ Lis n : Here Li was doped at the central Sn 
site and Vo was created at different O j sites, with 
.7 = 1-4 (see Fig. [Hb)) 

f) Vo+ Lioj: In this case, Vo was fixed and Li was 
doped at different O j sites, with j — 1,2,3 (see 
Fig. [He)) 

g) Vo+ Liint : Here Vo was fixed at the center of the 
supercell and Li was doped at interstitial sites. 



We performed calculations in the framework of den- 
sity functional theory (DFT) 31 using linear combi- 
nation of atomic orbitals (LCAO) as implemented in 
the SIESTA code — . The local density approximation 
(LDA)— was adopted for describing exchange-correlation 
interactions. We used standard norm-conserving pseu- 
dopotentials 34 in their fully nonlocal form — . Atomic 
positions and lattice parameters were optimized, us- 
ing a conjugate-gradient algorithm — , until the resid- 
ual Hellmann-Feynman forces converged to less than 0.05 
eV/A. A cutoff energy of 400 Ry for the real-space grid 
was adopted. 

The concentration of dopants and vacancies in a crystal 
depends upon its formation energies (Ef). The chemical 
potentials, which vary between the stoichiometric and 
diluted limits, depend on the material growth conditions 
and boundary conditions. We have calculated the forma- 
tion energies under equilibrium conditions, O-rich condi- 
tions and Sn-rich conditions. Under equilibrium condi- 
tions, for SnC>2 the chemical potentials of O and Sn sat- 
isfy the relationship /Lts n + 2/io = MSn0 2 i where /J.Sn0 2 j 
the chemical potential of bulk S11O2, is a constant value 
calculated as the total energy per Sn0 2 unit formula 
— . Under Sn rich conditions, ^sn=E(Sn motal ), fi — 
C 3 -" 2 ;^ ). O rich condition gives fi Q = (|)E(0 2 ), 

MSn = MSn02 — 2/Zo- 

The formation energies for systems with intrinsic de- 
fects, either Sn or O vacancies, can be calculated in the 
following way: 



E f = E d - E p + nfix, (1) 

where nx is chemical potential of X (X = Sn, O), n is the 
number of atoms removed from the system and E d and 
E p are the total energies of the defected and the pure 
system, respectively. 

The formation energy of the Li doped system can be 
calculated as 



E f = - (E d - E p + n/i x - m/XLi) , (2) 

n 

where fiu is the chemical potential of bulk lithium, cal- 
culated as total the energy per unit cell of bulk Li, n 
is number of atoms removed and m is number of atoms 
added to the system. 



III. RESULTS AND DISCUSSIONS 

In this section, we will discuss our DFT results. We 
will focus on Li-doped Sn02 systems, and then we will 
move to discuss the interactions of Sn and O vacancies 
with Li. 
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TABLE I: Formation energies calculated under equilibrium 
(E aq ), Sn-rich (i?sn) and O-rich (-Bo) conditions, in units 
(eV). The last column lists magnetic moments (MM) per su- 
percell , calculated in units of (j,b- Values in brackets show 
formation energies of the unrelaxed structures. 

System E cq E Sn E Q MM 

Li Sn -1.03(-0.89) 6.53(6.67) -1.03(-0.89) 3.00 

Lio -0.85 ( 5.72) -0.85 (5.72) 2.93(9.50) 0.00 

Li int -2.58(-1.63) -2.58(-1.63) -2.58(-1.64) 0.00 

V S n 6.87 (7.32) 14.44(14.88) 6.87(7.32) 4.00 

Vo 1.64 ( 1.69) 1.64 (1.69) 5.42(5.47) 0.00 



A. Li-doped SnCh 

First, we calculate, using the above equations, the de- 
fect formation energies of the doped SnC>2 systems, where 
Li was doped at the Sn, O, and interstitial sites, as stated 
before. The calculations (see Table |I])) show that the 
doping of Li at the Sn site has the lowest formation en- 
ergy (-1.03 eV) under stoichiometric and O-rich condi- 
tions. For comparison purposes, the formation energies 
of unrelaxed systems are also given in Table (JTJ) . It is clear 
that the structural relaxation is not very effective when 
Li is doped at the Sn site. This behavior can easily be 
understood in terms of the atomic sizes of the Sn and 
Li atoms. The atomic radii of Sn and Li are 1.41 Aand 
1.45 A, respectively. The doped Li is surrounded by six 
O atoms and the optimized bond lengths between them 
(LiO distances) have two different values, 1.97 A, and 
2.03 A, which are comparable to the pure bond lengths 
of SnO in Sn0 2 , i.e., 2.04 A and 2.08A. Since both Li 
and Sn approximately have the same atomic radius, then 
Li at the Sn site will not distort the structure and it is 
expected that the Sn site will be favorable for Li. Indeed 
our calculations also show that Li at the Sn site has the 
lowest formation energy, which agrees with experimental 
observations^. On the other hand, Li at the O site has a 
larger formation energy under the equilibrium condition. 
One can also see that the formation energy of Li at the O 
site is very large when structural relaxation is not allowed 
but decreases significantly when structural relaxation is 
permitted, which is followed by a large structural distor- 
tion. Our structural analysis shows that when Li is doped 
at the O site, it goes to interstitial site and leaves behind 
an O vacancy. So, it indicates that the formation energy 
of O is lowered mainly due to a structural deformation 
which is accompanied by the movement of Li to an in- 
terstitial position. At the same time, the system remains 
nonmagnetic. To be more confident, we also carried out 
separate calculations on Li at interstitial sites with and 
without O vacancies, and we got the same conclusion, 
which will be discussed below. Indeed, interstitial Li has 
the lowest formation energy, but it does not induce mag- 
netism in Sn02, so we will not focus on it. 

As stated above, doped Li prefers the Sn site under 




FIG. 2: (color online)The calculated total and partial density 
of states (DOS) of (a) Lis n , (b) Lio, and (c) Lii n t systems. 
Solid (red) and dashed (blue, green) lines represent s and p, 
d states, respectively. In the bottom panels the solid line 
represents the total DOS. The positive (negative) DOS shows 
majority (minority) spin states. The Fermi level (Ef) is set 
to zero. 



stoichiometric and O-rich conditions and has large mag- 
netic moment (3.00 /j,b). Since the nominal valence of Sn 
in a perfect Sn02 is Sn 4+ , and Li is a cation with valence 
of Li 1+ f22, when Li is doped at Sn site in pure Sn02, it 
donates one electron to compensate one hole among the 
four holes generated by the Sn deficiency. The three un- 
compensated holes, localized at the O sites, give a mag- 
netic moment of 3.00 //b per supercell. To see the atomic 
origin of magnetism in Li-doped Sn02 , we calculated the 
total and atom projected partial density of states(PDOS) 
of the Sn, O, and Li atoms (see Fig. The substitu- 
tional Li impurity polarizes the host band which gives 
rise to a induced impurity peak close to the valance band 
edge in the spin-down part of the total DOS, while the 
spin-up part is influenced slightly. Such a behavior sug- 
gests that Li behaves as a p type dopant in Sn02^ The 
low lying s orbitals of Li are spin-polarized and strongly 
hybridized with the p orbitals of O. The Fermi energy is 
mainly dominated by the p orbitals of O, which indicates 
that magnetism is mainly induced by the p orbitals and 
localized at the O atom. Indeed, the oxygen atoms sur- 
rounding the doped Li are the main ones that contribute 
to magnetism. We also see that the spin down state 
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is partially occupied which contributes to the magnetic 
moment. The majority s spin states of Li are completely 
occupied and the minority spin states are partially occu- 
pied, leading to a significant spin splitting near the Fermi 
level. 




FIG. 3: (color online)Spin density contours for the Li-doped 
Sn02 (HO) plane. Labels show atomic sites 



To further emphasize on the nature of magnetism, the 
spin density contours are shown in Fig. [3] It is interesting 
to see that Li polarizes the O and Sn atoms in opposite 
directions, and the polarization is not only limited to 
the nearest O atoms surrounding the Li atoms, but it 
also spreads to other O atoms which are far from Li. 
It seems that the origin of magnetism in light elements 
doped-Sn02 or Sn vacancies in Sn02 is the same i.e., 
the surrounding O atoms become polarized. However, 
the remarkable feature of Li doped SnC>2 is that there 
is a very small (negligible) induced magnetic moment at 
the Li site. This behavior is quite different from other 
light elements doped SnC>2 systems^— which suggests 
that Li behaves as a spin polarizer in SnC>2. We note 
that Li doped at O and interstitial sites does not induce 
magnetism (see Fig. Test calculations on a supercell 
of 2 x 2 x 3 (containing 72 atoms) were also carried out 
and we found the same behavior, i.e. Li at the Sn site 
gives 3.00 hb magnetic moment per supercell, whereas Li 
at the O site does not show any magnetic moment. 



B. Interaction of Sn vacancies with Li 

We will move now to discuss the interactions of intrin- 
sic defects (Sn and O vacancies) with Li. We found that 
a single Sn vacancy induces a very large magnetic mo- 
ment ~4.00 fiB which agrees with previous calculations 7 . 
However, Vs n has a very large formation energy in both 
the O-rich (6.87 eV) and Sn-rich(14.44 eV) conditions. 
Structural relaxation does not have a significant effect on 
the formation energy. To reduce the formation energy of 
Vsn, which is the core of this article, we considered an 
Sn vacancy at the center of the supercell, and doped Li 



TABLE II: Formation energies (in units of eV) of systems 
Vsn + Lisn, Vsn + Lio , and Vs n + Lii nt calculated under 
equilibrium (E eq ), Sn-rich (-Esn) and O-rich (Eo) conditions. 
The distance from Vs n to Lis n , Lio, and Lii nt is r — rj (in units 
of A) . The last column lists the calculated magnetic moments 
(MM) per supercell (in units of /ib). Values in brackets show 
the formation energies of the unrelaxed structures. 



System 


r - Tj 




771 


Eo 


MM 




V Sn 




f 01-7/1-7 oo\ 

6.87(7.32) 


1 A A A { 1 A C1Ci\ 

14.44(14.88) 


6.87(7.32) 4.00 


v Sn 


+ Li Snj 












V Sn 


+ Lisnl 


3.23 


1.55(2.19) 


9.11(10.39) 


1.55(2.83) 


7.00 


Vsn 


+ Li S n2 


"\ 7fi 


2.13(2.91) 


9.69(10.47) 


2.13(2.91) 


7.00 


V Sn 


+ Li Sn3 


4.80 


1.54(3.04) 


9.10(10.60) 


1.54(3.04) 


7.00 


VSn 


+ Lisn4 


5.78 


1.56(3.02) 


9.13(10.59) 


1.56(3.02) 


7.00 


VSn 


+ LisnS 


6.79 


2.80(3.04) 


10.36(10.61) 


2.80(3.04) 


7.00 


VSn 


+ Lisn6 


7.52 


2.91(3.14) 


10.48(10.70) 


2.91(3.14) 


7.00 


VSn 


+ Lioj 












V Sn 


+ Li i 


2.04 


-4.16(4.87) 


3.40(12.44) 


-0.38(8.65) 


1.00 


VSn 


+ Lion 


2.08 


-3.74(5.39) 


3.82(12.95) 


0.04(9.17) 


1.00 


VSn 


+ Lioin 


4.71 


6.09(10.91) 13.65(18.47) 


4.93(14.69) 1.00 


V S „ 


+ Li int 












V S „ 


+ Liint 


5.73 


1.83(4.04) 


9.40(11.60) 


1.83(4.04) 


3.00 



at different Sn sites (see Figure [TJ) . The formation en- 
ergies are given in Table (|TTJ) . By doping Li at the Sn 
site, Li significantly reduces the defect formation energy 
of Vsn- The case Vs n +Lis n 3 has the lowest formation en- 
ergy ~1.54 eV, which shows that Li doping is much more 
favorable energetically in SnC>2 than previously reported 
Zrj21 and Cr— doped-Sn02: The formation energies of 
the system with Zns n +Vs n are 16.50 eV and 7.00 eV 
under O-poor and the O-rich conditions^ 7 -, repectively. 
Doping Cr in SnC>2 with Vs n produces formation energies 
^12.00 eV and 2.00 eV under O-poor and O-rich condi- 
tions^, respectively. While in our case they are 9.10 eV 
and 1.54 eV under O-poor and O-rich conditions, respec- 
tively. It is interesting to see that Li not only stabilizes 
Vsn, but also has the largest magnetic moment 7.00 //b 
per supercell. When the Li atom occupies the Sn site, 
it donates one electron to compensate one hole among 
the four holes generated by the substituted Sn. These 
three uncompensated holes along with the four holes of 
neutral Vs n give a total of 7.00 /ib magnetic moment per 
supercell. 

Once found that Li at Sn sites significantly reduces the 
defect formation energy of Sn vacancies, then we doped 
Li at the O site in the presence of Vs n - In this case, we 
also fixed Vs n and doped Li at different O sites, mainly 
considering those O atoms which are near and far from 
Vsn- The calculated formation energies are shown in Ta- 
ble ((TTJ) . It is clear that the unrelaxed systems have large 
positive formation energies (endothermic processes) , but 
after structural relaxation we found negative formation 
energies (exothermic processes). When Li was doped at 
O sites around Vs n it goes to an interstitial site (near 
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^ Vv#H ^ 



FIG. 4: (color online) Schematic representation of structural 
relaxations in the Lio + Vs n system. Green, yellow, and red 
balls represent Sn, O, and Li atoms, respectively. The Sn 
vacancy (Vsn) is represented by a green dashed circle. The 
arrow represents movement of the Li atom towards Vs n during 
structural relaxation. 



Vsn) to form a stable structure, after structural relax- 
ation. This process leaves behind an oxygen vacancy and 
Li interstitial near Vs n (see Fig.[4j). This exothermic pro- 
cess further indicates that doped Li does not prefer the 
O site, rather it prefers the interstitial site, as stated 
above. As shown also above, the Lio and Lii nt systems 
have no magnetic moment, whereas Vs n has large mag- 
netic moment. However, when Li is doped at the O site 
in a system with Vs n , it breaks the bonds with the neigh- 
boring Sn atoms and goes to the nearest interstitial site 
leaving behind an O vacancy. Three out of the four holes 
generated by Vs n are compensated by 2~ charge state of 
Vo and one electron of Li, leaving behind a single hole 
that gives a magnetic moment of 1.00 /jb- When Li is far 
from Vsn (Vsn+Lio/// case), the Li atom does not form 
a stable structure. 

We discuss now the interaction of Sn vacancies with 
Li; nt . In this case, we doped Li at an interstitial site 
which was 5.37 A away from Vs n - The defect formation 
energy of this system is given in Table |Tj]) . Although the 
Liint + Vs n system behaves endothermically, the forma- 
tion energy of Vs n decreases so much as compared with 
Vsn system without Li; nt . Note that the structural re- 
laxation does not change the sign of Ef . The Lij nt + Vs n 
system also shows a large magnetic moment, 3.00 /ib per 
supercell, or in other words, when Li is at the interstitial 
site it interacts with Vs n and reduces the magnetic mo- 
ment of Vs n . From the above thermodynamics of differ- 
ent defects in SnC>2 , we can conclude that among all these 
defects the Vs n + Lio complex has the lowest energy. In 
the presence of vacancies (Sn and O), Li prefers to oc- 
cupy the O site rather than going into Sn or interstitial 
sites. Although, interstitial Li also reduces the formation 



energy of Vs n , it is not a favorable process. Note that 
different defective systems of Li-doped ZnO showed the 
same type of behavior in a previously reported work<2£ 
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FIG. 5: (color online)The calculated total and partial den- 
sity of states (DOS) of (a) Vsn+Lisn, (b) Vsn+ Lio, and (c) 
Vsn+Liint systems. Solid (red) and dashed (blue, green) lines 
represent s and p, d states, respectively. The bottom panel 
solid line represents the total DOS. The positive (negative) 
DOS shows majority (minority) spin states. The Fermi level 
(Ef) is set to zero. 



It is very interesting that irrespective of the location 
of Li in Sn defective systems, it always reduces the for- 
mation energy of Vs n and the systems remain magnetic. 
To further understand the electronic structure of SnC>2 
and the interaction of Li with native defects, we show the 
calculated total and atom projected densities of states in 
Fig. [5] The interaction of Li with native defects mainly 
affect the local DOS of the O atoms. Also, Li-0 type of 
hybridizations can be seen around the Fermi levels. In all 
cases, the oxygen spin-up p orbitals are occupied whereas 
the spin-down p orbitals are partially occupied, similar 
to the DOS of V Sn - 7 It is noticeable that all the DOS fig- 
ures show half metallic behavior except Vsn + Lio which 
anyhow has nearly half metallic behavior, mainly due to 
the p orbitals of the O atoms. It is also noticeable that 
when Li goes to an interstitial site in the presence of a Sn 
vacancy, the oxygen p orbitals are no longer degenerated 
at the Fermi level. 
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Interaction of O vacancies with Li 



TABLE III: Formation energies (in units of eV) of systems 
Vo + Lisn , Vo + Lio , and Vo + Liint calculated under equi- 
librium (Eoq), Sn-rich (Esn) and O-rich (-Eo) conditions. The 
distance from Vo to Lig n , Lio, and Lii nt is r — r, (in units of 
A). The last column lists the calculated magnetic moments 
(MM) per supercell (in units of Mb)- Values in brackets show 
the formation energies of the unrelaxed structures. 



System 


r-rj 


-Ecq 


Esn 


E MM 


Vo 




1.64( 1.69) 


1.64( 1.69) 


5.42(5.47) 0.00 


Voj + Li Sn 










Voi + Li Sn 


2.04 


-2.55(-0.07) 


5.01(7.49) 


1.23(3.71) 1.00 


V02 + Li Sn 


2.08 


-2.34(-0.50) 


5.22(7.06) 


1.44(3.28) 1.00 


V03 + Li Sn 


4.71 


-3.69(-2.11) 


3.87(-2.11) 


0.09(1.67) 1.00 


V04 + Lisn 


5.71 


-2.09(-2.69) 


3.39(4.86) 


-0.39(1.08) 1.00 


Vo + Lioj 










Vo + Lioi 


2.94 


0.15(3.48) 


0.15(3.48) 


3.93(7.27) 0.00 


Vo + Li 2 


3.47 


0.21(5.01) 


0.21(5.01) 


3.98(8.79) 0.00 


Vo + Lio3 


5.15 


0.30(3.80) 


0.30(3.80) 


4.08(7.59) 0.00 


Vo + Liint 










Vo + Liint 


3.48 


-0.96(0.37) 


-0.96(0.37) 


1.36(4.16) 0.00 



It is believed that Sn vacancies have higher formation 
energy as compared with O vacancies. Indeed, we found 
that the O vacancy has defect formation energy (1.64 eV) 
much smaller than that for the Sn vacancy. We therefore 
considered also the interactions of O vacancies with Li. 
First, we doped Li at the Sn site and varied the posi- 
tion of the O vacancy. The calculated defect formation 
energies are shown in Table (|IIIj) . Interestingly, doping 
of Li at the Sn site appreciably reduced the defect for- 
mation energy of Vo- The defect complex V03 + Lig n 
has the lowest formation energy, -3.69 eV in equilibrium 
conditions. Such a favorable process implies that oxygen 
vacancies also promote magnetism in doped systems. We 
must stress that oxygen vacancies by itself are nonmag- 
netic^, but however, become magnetic in the presence of 
Lisn, which is magnetic. Either Lig n or Vo have large 
defect formation energies, but the defect complex Vo + 
Lisn bas a lower energy than the individual defects. It 
is interesting to see that these defects complex always 
give magnetism with magnetic moments of 1.00 /ib, no 
matter how far is Vo from Lig n . Such behavior is differ- 
ent from Co-doped SnC>2 where it has been shown that 
Vo quenches magnetism when Cos n is away from Vo^ 
On the other hand, Ni-doped Sn02 does not show mag- 
netism without Vo 42 We therefore believe that Li-doped 
SnC>2 is another candidate material for spintronics. We 
also found that Vo and Lio do not induce magnetism in 
SnC-2, as seen in Table (|IIII) . but the formation energy 
of this complex is always bigger than that of the Vo + 
Lisn complex, excluding under Sn-rich conditions. There 
are some experimental reports which claim that oxygen 
vacancies are responsible for room temperature ferromag- 



netism in oxides^ but to date, there are no experiments 
reports which show that a single neutral Vo is the main 
source of magnetism in Sn02- Generally, it is believed 
that Vo promotes magnetism in doped Sn02 i 39 ' 40 which 
is proven by our extensive calculations (see Table (|IIip). 
We also considered Li at different O sites in the presence 
of O vacancies but the calculated E[ show that these con- 
figurations are less stable than the Vo + Lisn (excluding 
Sn-rich conditions). 

Finally, the interaction of O vacancies with Lij nt was 
also taken into account. The calculated E[ of Li; nt +Vo 
are given in Table pil[) which shows that the behavior of 
Liint +Vo is quite different from Lii„t+Vs n - Liint reduces 
the formation energy of Vs n , but however Lij nt does not 
decrease the formation energy of Vo- Also such kinds 
of defect complexes have zero magnetic moments. From 
Table (JTTJ) and Table (|III[) . we can see that Ef can be re- 
duced by Li, but at the same time thermodynamics also 
tell us that in the presence of a Sn vacancy the probable 
location of Li will be the O site where structural relax- 
ation is also expected. On the other hand, Li prefers the 
Sn site rather than the O site, when the interactions of 
O vacancies with Lio is considered. 
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FIG. 6: (color online)The calculated total and partial den- 
sity of states (DOS) of (a) Vo+Lis n (b) Vo+Lio, and (c) 
Vo+Liint systems. Solid (red) and dashed (blue, green) lines 
represent s and p, d states, respectively. The bottom panel 
solid line represents the total DOS. The positive (negative) 
DOS shows majority (minority) spin states. The Fermi level 
(Ef) is set to zero. 
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Before summarizing our work, the electronic structure 
of the above mentioned defect complexes will be pre- 
sented. FigurelUa) shows that the electronic structure is 
similar to Lis n , i-e., the magnetism is mainly contributed 
by the O p orbitals and the Li-0 type hybridization can 
also be seen near the Fermi level. No magnetism can 
be seen in Fig.^b) which represents the Vo+Lio defect 
system. The impurity peak around 2 eV below E-p is 
contributed by the Li atom. The conduction and valance 
bands are separated from each other and the Fermi level 
crosses the conduction band, which is mainly formed by 
Sn s and Li orbitals. This electronic structure indicates 
that Li in the Vo system will behave as an electron 
donor. The electronic structure of Vo+Li;„t [Fig. EJc)] 
also shows nonmagnetic behavior. The O p orbitals are 
completely occupied and are far from Ep. The Li and 
Sn atoms occupy the conduction band, i.e., Vo+Lii n t 
is also a n type material. Comparing the DOS of the 
Liint system with that of the Vo+Li; n t, we see that Vo 
pushes the valance band of Li int to lower energy states. 
From the given DOS in Figs. @, ©, and ([5]), we may 
say that those systems in which E-p cuts the conduction 
bands, i.e., which behave as electron donors, do not show 
magnetism. This indicates that magnetism in Sn02 is 
mediated by holes and destroyed by electrons. Further 
experimental work is needed to validate our theoretical 
calculations. 



at Sn, O, and interstitial sites and it was shown that it 
can induce magnetism in Sn02 when doped at the Sn site. 
No magnetism was found however when Li was doped at 
the O and interstitial sites. The defect formation ener- 
gies showed that Li doped at the Sn site is more favorable 
than doped at the O site. We found that Li at the Sn 
site also shows a large magnetic moment (3.00/ib), and 
the origin of magnetism was discussed in terms of elec- 
tronic structures and spin densities. Native defects, Sn 
vacancy (Vs n ) and O vacancy (Vo), were also studied 
and it was observed that magnetic Vs n has higher defect 
formation energy than non magnetic Vo- To reduce the 
defect formation energies of these native defects, the in- 
teractions of Li with Vsn and Vo were also studied. The 
calculated defect formation energies of these native de- 
fects were significantly decreased by Li. Our calculations 
showed that Vo helps to promote magnetism in Li doped 
Sn02, in agreement with general experimental specula- 
tions. Structural relaxations were shown to be important 
when the interactions of Li with Vs n were considered. 
Comparison of our studied system with the previously 
synthesized systems was also discussed and it was con- 
cluded that Li doped Sn02 is another good candidate in 
the filed of magnetism. 
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